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ABSTRACT 

High-resolution (R ss 100 000), high signal-to-noise spectra of M71 giants have been obtained with 
HIRES at the Keck I Telescope in order to measure their Mg isotopic ratios, as well as elemental 
abundances of C, N, O, Na, Mg, Al, Si, Ca, Ti, Ni, Zr and La. We demonstrate that M71 has two 
populations, the first having weak CN, normal O, Na, Mg, and Al, and a low ratio of 26 Mg/Mg 
(~4%) consistent with models of galactic chemical evolution with no contribution from AGB stars. 
The Galactic halo could have been formed from the dissolution of globular clusters prior to their inter- 
mediate mass stars reaching the AGB. The second population has enhanced Na and Al accompanied 
by lower O and by higher 26 Mg/Mg (~8%), consistent with models which do incorporate ejecta from 
AGB stars via normal stellar winds. All the M71 giants have identical [Fe/H], [Si/Fe], [Ca/Fe], [Ti/Fe] 
and [Ni/Fe] to within a = 0.04 dex (10%). We therefore infer that the timescale for formation of the 
first generation of stars we see today in this globular cluster must be sufficiently short to avoid a con- 
tribution from AGB stars, i.e. less than ~0.3 Gyr. Furthermore, the Mg isotopic ratios in the second 
M71 population, combined with their elemental abundances for the light elements, demonstrate that 
the difference must be the result of adding in the ejecta of intermediate mass AGB stars. Finally we 
suggest that the low amplitude of the abundance variations of the light elements within M71 is due 
to a combination of its low mass and its relatively high Fe-metallicity. 

Subject headings: stars: abundances - stars: atmospheres - stars: evolution - globular clusters: M71 



1. INTRODUCTION 

Galactic globular clusters (GCs) show star-to-star vari- 
ations in the light elements including Li, C, N, O, Na, 
Mg, and Al discovered more than 30 years ago. Corre- 
lations and anti-correlations among O, Na, Mg, and Al 
are seen in essentially all GCs, but not in field stars, as 
reviewed bv lGratton. Sneden fc Carrettal (|2004f h These 
correlations are highl y suggestive of proton burning of H 
at high temperatures (jDenisenkov fc Dcniscn kovalll990f ). 
However, detection of these correlations among GC main 
sequence stars which have not yet evolved sufficiently for 
their central temperatures to become hot enough for the 
necess ary nucle ar reactions to ope r ate b y iGratton et al 
(|2001h and by iRarmrez fc Cohenl (|2002ft demonstrated 
definitively that these patterns are imprinted from an 
external source; they cannot be the result of internal nu- 
cleosynthesis follo wed by mixing to the su r face o f pro- 
cessed material. ICohen. Brilev fc Stetson! ((20051 ) have 
shown that superposed on this is the imprint of mix- 
ing of material processed through H burning via the CN 
cycle, but only among the most luminous of the globular 
cluster giants. 

The scenario generally invoked to explain this is that 
more than one generation of stars was formed in the 
GC; first one that produced the heavy elements we see 
today within a system where the gas was well mixed 
throughout the entire volume. This was followed by 
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Technology, Pasadena, California 91125 (jlc@astro.caltech.edu) 



the second generation of stars, including low mass stars 
that survive to the present, then a subsequent stellar 
generation formed from gas contaminated by the ejecta 
of an earlier generation that was not uniformly mixed 
throughout the volume. It is generally presumed that 
material processed within the interiors of AGB stars, 
mixed to the stellar surface, then ejected by stellar winds, 
is the culprit producing the star-t o-star variations de - 
tected among GC stars, s ee, e .g. iFenner et all (l2003tl , 
ICohen. Brilev fc Stetson] (|2005f L I Ventura fc DAntonal 
(2008), and references therein. Recently rapidly rotat- 
ing massive stars have also been invoked for this purpose 
as during their main sequence phase they transport ma- 
terial processed through H burning at high t empe ratures 
to their surfaces, see, e.g. iDecressin et a 1 (l2007h . This 
material has an escape velocity low enough to be confined 
to the cluster, unlike that of the subsequent radiatively 
driven winds or the eventual SN, where the ejected mate- 
rial has such a high velocity that it is lost by the cluster. 

However, some of the details of the behavior seen thus 
far do not match the predictions for high temperature 
H-burning. The most discrepant case is that of the Mg 
isotopes. The first atte mpt to measure Mg isotopic ratios 
for GC giants wa s bv IShetrqnd (Il996f) . wit h additional 
recent studies bv lYong et al.1 (|2003f ) and bv lYong et al.1 
( 2006). All of these authors found substantial contribu- 
tions from the heavier Mg isotopes, reaching up to half 
of the total Mg, far larger than the solar ratio, which 
is 79:10:11 for 24 Mg: 25 Mg: 26 Mg. Measurement of the 
Mg isotopic ratios requires detection of weak contribu- 
tions from the rarer heavier Mg isotopes in the wings of 
much stronger 24 MgH lines, hence exquisite spectra of 
very high spectral resolution (^100,000) and very high 
signal-to-noise ratio, as well as careful attention to de- 
tails of the line list used for the spectral synthesis and to 
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the analysis procedure. 

Here we present a new attempt to measure the Mg iso- 
topic ratios in the giants of the relatively metal-rich GC 
M71 which we hope fulfills the criteria described above to 
achieve a reliable result. We exploit our exquisite spec- 
tra to derive precision abundances for several additional 
elements in M71, including the light elements O, Mg, 
Na, and Al. We discuss our results in the context of the 
history of star formation in globular clusters, and the 
contribution of AGB stars to GCs and to the halo field. 

2. OBSERVATIONS 

The sample consists of 9 giant stars in the globular 
cluster M71, eight observed with HIRES (Vogt et al. 
1994) at the Keck I telescope by us and one (M71 A4) 
observed with HDS at the Subaru telescope by Yong et 
al. (2006) 4 . 

The Keck observations were obtained in August 2004, 
June 2005 and September 2007. All our data was ob- 
tained after the HIRES upgrade (August 2004), thus 
taking advantage of improvements in efficiency, spectral 
coverage and spectral resolution. A resolving power of 
B w 10 5 was achieved using a 0.4"-wide slit, accepting 
a substantial light loss at the slit in order to achieve the 
necessary spectral resolution. The signal-to-noise level 
exceeded 100 per pixel, with 1.3 km s _1 /pixel. In Fig. 1 
we compare our HIRES spectrum of M71 1-77 with the 
HDS spectrum of M71 A4 (Yong et al. 2006), which is 
of similar quality (R w 90 000). 

The spectral orders were extracted using both IRAF 
and MAKEE 5 . Further data reductions (Doppler correc- 
tion, continuum normalization, and combining spectra) 
were performed with IRAF. 

3. ANALYSIS 

Initially stellar parameters were obtained from 
photometry, followin g the precepts set out in 
ICohen. Behr fc Brilevl (|2001h . We are fortunate 
to now have J and K s photometry from 2MASS 
(jSkrutskie et all 119971 : ICutri et all 120031 ) (not available 
at the time of the initial M71 abundance study) and 
optical photometry with mo r e exte nsive coverage of the 
area of M71 from I Stetson! (|2000f ) as updated on his 
on-line database available through the CADC. We adopt 
E(B - V) = 0.25 mag from the Feb. 2003 version of the 
online database of Harris (1996) 6 . 

Our preliminary analysis using the photometric stellar 
parameters and Kurucz overshooting model atmospheres 
(Castelli, Gratton & Kurucz 1997) revealed that M71 
seems to show only very slight abundance variations (if 
any) in some elements (e.g. Mg), and therefore very high 
precision is needed to ascertain if abundance variations 
are present or not. Seeking the highest possible precision, 
and given the potential presence of reddening variations 
across the field of M71, we decided to perform the anal- 
ysis using spectroscopically determined stellar parame- 
ters. The stars have similar temperatures (T e gp = 4045 

4 The spectrum of M71 A4 has been kindly made available to us 
by D. Yong and W. Aoki 

5 MAKEE was developed by T. A. Barlow specifically 
for reduction of Keck HIRES data. It is freel y availabl e at 
http: / / www2.keck.hawaii.edu/realpublic/inst/hires / data_reduction. 

6 http:/ /physwww. monaster. ca/~harris/Databascs. html 



± 225 K), so this approach should minimize errors in the 
relative stellar parameters. 

Excitation and ionization balances of Fei and Fell 
were performed with the 2002 version of MOOG (Sneden 
1973), using the Fei line list of Alves-Brito et al. (2009, 
in preparation) and the Fe II line list of Hekker & Me- 
lendez (2007), which is described in Melendez & Barbuy 
(2009). These line lists have been carefully scrutinized 
to avoid significant blends (e.g. CN) in metal-rich cool 
giants. All equivalent widths were measured by hand us- 
ing IRAF, employing the deblend option when necessary. 
Microturbulence was obtained by flattening the trend in 
the iron abundance from Fe I lines versus reduced equiv- 
alent width. The zero-points of our spectroscopic stel- 
lar parameters have been determined by Melendez et al. 
(2008) using many bright field giants with precise effec- 
tive temperatures using the infrared flux method tem- 
perature scale of Ramirez & Melendez (2005), and with 
surface gravities determined with Hipparcos parallaxes 
and Yonsei-Yale (Demarque et al. 2004) and Padova 
isochrones (da Silva et al. 2006). 

Once the stellar parameters (T e g , log g, v m i C ) were set, 
the macroturbulence (u mac ) was determined employing 
five Fei lines (605.60, 607.85, 609.67, 612.02, 615.16 nm) 
and checked with other lines present in the same region 
(605 to 615 nm). Note that in previous works on Mg iso- 
topic ratios (e.g. McWilliam & Lambert 1988; Yong et al. 
2003) only the Ni I 511.54 nm and the Ti I 514.55 nm were 
used to obtain the macroturbulent broadening. However, 
we find that these lines arc not useful in metal-rich cool 
giants, because the region around each of them is badly 
blended. The region we have chosen to use for this pur- 
pose (605-615 nm) is much cleaner, although we have also 
included in our spectrum synthesis blends by atomic, C2, 
and CN lines. We recommend the use of this region to 
determine v m ac in future studies of Mg isotopic ratios in 
relatively metal-rich very cool stars. 

The line lists of Mg H and C2 (0-0 band ) have been 
previously described in IMelendez fc Cohe"nl ([2007) , and 
the exten sion of our C2 line li s t to o ther bands is de- 
scribed in lMelendez fc Asplundl (12008ft. Th e CN l ine list 
adopted here is from IMelendez fc Barbuvl |T999). For 
the atomic lines, we have used when available transition 
probabilities based on laboratory measurements, other- 
wise we have adopted astrophysical ^/-values based on 
Arcturus. We only used the best, most isolated atomic 
lines within the wavelength range covered by our spectra, 
listed in Table Q] 

The Mg isotopic ratio s were determined b y spe ctral 
synthesis, as described in lMelendez fc Cohenl (|2007h . ex- 
cept that in cool metal-rich giants the MgH feature at 
513.46 nm (usually employed for Mg isotopic ratios; e.g. 
Gay & Lambert 2000) seems blended, so instead we use 
the 513.43 nm feature. Two other MgH features at 513.87 
and 514.02 nm (e.g. Yong et al. 2003; Melendez & 
Cohen 2007) were also used for obtaining the Mg iso- 
topic ratios. After the first trials, it was clear that the 
24 Mg: 25 Mg: 26 Mg ratios were no larger than terrestrial 
(79:10:11), therefore we have computed synthetic spec- 
tra with isotopic ratios ranging from 24 Mg: 25 Mg: 26 Mg = 
100:0:0 to 76:12:12. 
The isotopic ratios were determined by \ 2 fits, by com- 
S(Oj — Si) /a 2 , where Oi and Si represents 
he observed and synthetic spectrum, respectively, and 
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a = (S/N) -1 . Examples of the x 2 fits are shown in 
Figs. 2 and 3 for the three recommended MgH features. 
The resulting Mg isotope ratios are presented in Tab. 2. 
In this table are also given the standard deviation from 
the three different MgH features, which for 25 Mg and 
26 Mg have a typical value of a = 1.7% and 1.9%, re- 
spectively, corresponding to standard errors of 1.0% and 
1.1%. As discussed below, the true error for 26 Mg is actu- 
ally somewhat smaller than for 25 Mg. Due to the larger 
isotopic separation of the 26 MgH lines, the determination 
of 26 Mg/ 24 Mg is more reliable than 25 Mg/ 24 Mg. 

The elemental abundances of C, N, O, Na, Mg and Al 
were obtained by spectral synthesis, taking into account 
blends by atomic, C2 and CN lines, and the abundances 
of Si, Ca, Ti, Ni and Zr were obtained from equivalent 
widths. Carbon abundances were obtained from C2 lines 
around 563 nm and checked using CH lines around 430 
nm (Plcz ct al. 2008, see also Plez & Cohen 2005). Ni- 
trogen abundances were obtained using CN lines around 
630-637 nm and 669-671 nm. The CN-rich stars have 
N abundances ~3 times higher than the CN-weak stars, 
but the difference may be higher because TiO blends pre- 
vent a precise determination of N abundances in the CN- 
weak stars. The elemental abundance ratios are given in 
Tabs. 2 and 3. 

The iron abundances, [C,N,0/Fe] abundance ratios, 
and the C+N and C+N+O abundance sums, are shown 
in Fig. 4 as a function of effective temperature. A simi- 
lar plot for the isotopic 24 < 25 < 26 Mg/Mg, 25 < 26 Mg/ 24 Mg ra- 
tios and the elemental abundance ratios [Na,Mg,Al/Fc] 
is presented in Fig. 5, and the [Si,Ca,Ti,Ni,Zr,La/Fe] 
ratios are shown in Fig. 6. 

In these and subsequent figures CN-weak and CN-rich 
giants are represented by open and filled circles, respec- 
tively. The CN status of our sample stars have been 
obtained from Smith & Norris (1982) and Lee (2005), 
except for star M71 I, for which no previous information 
on its CN bands is available in the literature. Since this 
star shows strong CN bands in our HIRES spectrum, we 
have also included it as a CN-strong star. 

As discussed below, some abundance ratios show 
trends with T ff . The coefficients of the fits for the CN- 
weak stars are given in Tab. 4. 

4. DISCUSSION 

4.1. Correlations between 24 < 25 < 26 Mg/Mg, 0, Na, Mg 

and Al 

The Mg isotopic ratios and the elemental abundances 
of O, Na, Mg and Al for the CN-weak stars seem to 
show some trend with effective temperature, probably 
due to NLTE and 3D effects (Asplund 2005; Collet et al. 
2007). After correcting for trends with T e ff, the scatter 
in [0,Na,Mg,Al/Fe] for the CN-weak stars is only 0.018, 
0.039, 0.018 and 0.021 dex, respectively. The CN-rich 
stars depart from the behavior of the CN-weak stars, 
showing larger 25 ' 26 Mg/Mg ratios, larger Na and Al ele- 
mental abundances, and lower O and Mg elemental abun- 
dances. 

In Figs. 7-11 are shown the 24 Mg/Mg, 25 Mg/Mg, 
25 Mg/ 24 Mg, 26 Mg/Mg and 26 Mg/ 24 Mg ratios as a func- 
tion of [0,Na,Mg,Al/Fe]. In these plots the trend with 
effective temperature (Figs. 4-5, Tab. 4) has been cor- 
rected in the [0,Na,Mg,Al/Fe] ratios, to prevent spurious 
trends. The 25,26 Mg/ 24 Mg isotopic ratios show a larger 



spread than 25 ' 26 Mg/Mg due to the spread in 24 Mg/Mg 
(CN-rich giants show lower 24 Mg/Mg than CN-weak gi- 
ants). 

In Fig. 7 we show that 24 Mg in CN-rich stars is an- 
ticorrelated with Na and probably also with Al. The 
25 Mg/Mg ratios are less reliable than the 26 Mg/Mg ra- 
tios. Although Table 2 suggest higher uncertainties for 
26 Mg than for 25 Mg, with standard errors of 1.1% and 
1.0%, respectively, the star-to-star scatter for the CN- 
weak stars actually shows that the errors for 25 Mg and 
26 Mg are 0.9% and 0.8%, respectively, i.e., in both cases 
somewhat lower than the errors estimated from Table 2. 
The higher error bar (±0.9%) for 25 Mg is expected due 
to the smaller isotopic shift of the 25 MgH lines. Hence, 
25 Mg does not show clear trends besides the fact that 
25 Mg is enhanced in CN-rich giants (Figs. 8-9). The 
CN-weak stars do not show any correlation between 26 Mg 
and the elemental abundances of Na, O, Mg and Al (Figs. 
10-11), but the CN-rich giants show strong, weak and no 
correlations between 26 Mg and the elemental abundances 
of Na and Al, O, and Mg, respectively. 

Fig. 12 shows the correlations between the elemental 
abundances of O, Na, Mg and Al, which have been cor- 
rected for trends with T e g (Figs. 4-5, Tab. 4), and there- 
fore any remaining trend or scatter should be real. Inter- 
estingly, even though the CN-weak stars show only small 
scatter in their abundances ratios (0.02-0.04 dex), and 
seem to be roughly constant between the uncertainties, 
there is a hint of a correlation between Na and Al. On 
the other hand, Al and Mg, and O and Na, seem to be 
anti-correlated, although the evidence is weaker for the 
0:Na anti-correlation. The CN-strong stars show strong 
correlation and anti-correlation between Na and Al, and 
O and Na, respectively. 

4.2. C+N and C+N+O in M71 giants 

As can be seen in Fig. 4, the C+N abundance sum 
is larger in CN-rich than in CN-weak M71 giants. This 
is mainly due to the large N enhancement of the CN- 
rich giants, which have N abundances ~0.5 dex higher 
than the CN-weak giants. Keck low resolution spectra 
of CN-weak and CN-strong main sequence stars in M71 
show that CN-strong dwarfs are also enhanced in nitro- 
gen when compared to CN-weak dwarfs (Briley & Cohen 
2001). 

The C+N+O abundance sum is constant within 0.1 
dex (Fig. 4), in agreement with other high resolu- 
tion analysis of GC giants in the literature, which find 
C+N+O constant within 0.3 dex for NGC 6712 ([Fe/H] 
= -1.0, Yong et al. 2008), M4 ([Fe/H] = -1.1, Smith et 
al. 2005) and M13 ([Fe/H] = -1.5, Cohen & Melendez 
2005). 

Recently, Yong et al. (2009) have found a large 
C+N+O spread of 0.57 dex for NGC 1851 giants ([Fe/H] 
= -1.2). This spread exceeds their estimated uncertainty 
(0.14 dex), and is much larger than that observed in other 
GCs, where no significant spread is found. Yong et al. 
(2009) interpret the large spread in C+N+O as a signa- 
ture of AGB pollution in NGC 1851, and they conclude 
that if the AGB scenario is applicable to other GCs with 
constant C+N+O, then perhaps the masses of the AGB 
polluting stars in NGC 1851 were lower than in other 
GCs. This is probably also the reason why NGC 1851, 
unlike other GCs, shows a large spread in the s-process 
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elements Zr and La (Yong et al. 2009). 

4.3. Chemical evolution of the Mg isotopes 

As discussed in lMelendez fe Cohen] (|2007l ). the lightest 
and most abun dant Mg isotope can be fo rmed in mas- 
sive stars (e.g. IWooslev fe Weaver] [19 95;) as a primary 
isotope from H. The heavier isotopes are formed as sec- 
ondar y isotopes, as well as in int ermediate-mass AGB 
stars (Karaka s fe Lattanziol 120071 ). so the isotopic ra- 
tios 25 ' 26 Mg/ 24 Mg increase with the onset of AGB stars. 
Therefore Mg isotopic ratios in halo stars can be used as 
a chronometer to constrain the rise of the AGB in our 
Galaxy. In our earlier paper we applied this technique 
to a sample of metal-poor Galactic halo field dwarfs to 
establish that the AGB did not make a significant con- 
tribution to Mg in the Galactic halo at least until [Fe/H] 
~ — 1.3 dex was reached. 

Fig. 13 shows the chemical evolution of 26 Mg/ 24 Mg 
for halo field dwarfs (Melendez & C ohen 2007), adding 
in the M71 giants, with models from Fcn ner et all (|2003f ) 
for the solar neighborhood both with and without the 
contribution of AGB stars. We see that the M71 CN- 
weak giants have very low Mg isotopic ratios 26 Mg/Mg 
(~4%) consistent with models of galactic chemical evolu- 
tion with no contribution from AGB stars, and extending 
our previous results based on metal-poor halo field dwarfs 
to still higher metallicity. Base d on the calculations 
of ste llar yields for AGB stars b y Karakas fe Lattanziol 
(2003), updated in lKarakas fe Lattanziol (|2007l ). stars of 
initial mass 3-6 Mq are required to contribute signifi- 
cant amounts of the heavier Mg isotopes. Stars in this 
mass range reach the upper AGB at an age of ^0.3 Gyr; 
see, e.g., the Padova evolutionary tracks 7 . For a GC 
as metal-rich as is M71, this timescale must reflect the 
minimum age difference between the CN-weak and the 
CN-rich stars in the GC as well as the maximum range in 
age of the CN-weak stars we see today in M71, assuming 
uniform mixing of the gas within the GC. Furthermore, 
the enrichment of the halo gas from which M71 eventu- 
ally formed must have reached such a high Fe-metallicity 
without a substantial contribution from proton burning 
of H at high temperatures. We argue that AGB stars 
are the source of the "polluted" material, rather than it 
coming from rapidly rotating massive stars, on the ba- 
sis of timescales, which would be uncomfortably short in 
the latter case, and also the uniformity of the heavy ele- 
ments within all stars of M71. Our precision abundances 
define an upper limit to the range of abundances of the 
heavy elements in M71; all the M71 giants have identical 
[Fe/H], [Si/Fe], [Ca/Fe], [Ti/Fe] and [Ni/Fe] to within 
a = 0.04 dex (10%) after the trends with T e ff are re- 
moved (see Figs. 4 and 6). This would be hard to achieve 
if massive stars were the source of polluting gas in GCs. 

On the other hand, the M71 giants with strong CN (ac- 
companied by lower O with higher Na and Al than those 
of the weak CN giants) show higher Mg isotopic ratios 
26 Mg/Mg (^8%), but these are st i ll quit e low compared 
to previous studies by lYong et all (|2006l ) mostly in GCs 
of lower metallic ity. In parti c ular, star M71 A4 has been 
also analyzed bv lYong et alJ (|2006h . who find higher iso- 
topic ratios than ours. This must be due to the different 
analysis techniques, as the same observed spectrum was 

7 http://pleiadi.pd.astro.it/ 



used. In particular, note that their macroturbulent ve- 
locity for this star is lower. As shown in Fig. 13, our 
newly measured Mg isotope ratios for the CN-rich giants 
in M71 lie on the predicted curve for models which in- 
clude the appropriate contribution to the chemical inven- 
tory from material processed through high temperature 
H-burning which we believe comes from AGB stars dur- 
ing the normal course of evolution of intermediate mass 
stars. 

4.4. Neutron capture processes 
Fig. 6 displays [La/Fe] vs T e g for the M71 sample of 
luminous giants. The 6390 A La II line was synthesized 
with hyperfine stru cture and transition probability from 
lLawler et al] |2001) . The results are more uncertain than 
for most of the other elements discussed here as there is 
only one La II line which was sufficiently unblended to 
be used, and it lies within a CN band. However, there 
appears to be a slight excess of La among the CN-strong 
M71 giants. This excess is small enough that it would 
n ot have been detect e d in p revious analyses such as that 
of lRamirez fe Cohen! (|2002l ). whose spectra were of lower 
spectral resolution and SNR, with consequent larger un- 
certainties. 

La, together with Ba 8 , are the archetypical elements 
indicating a contribution to the star's chemical inven- 
tory via s-process neutron capture. But this provides no 
discrimination in the source of the s-process, be it in- 
termediate mass AGB stars or young massive rotating 
stars. The former, especially in metal-poor AGB stars, 
is we ll known as a site for such production (|Busso et all 
l200lh. while very recent work by the NuGrid project 
( Hirschi et al.ll2008l ) demonstrates that the s-process can 
operate efficiently in rotating young massive stars. 

Zr is also believed to be produced primarily in AGB 
stars via the s-process of neutron capture on Fe seeds. 
However, we do not see any detectable difference for the 
mean [Zr/Fe] between the CN-rich and the CN-weak M71 
giants, while differences in the isotopic ratio 26 Mg/Mg 
are apparent, and we believe there are differences in 
[La/Fe] as well. There are four detected Zr I lines, so the 
lack of detectable star-to-star variation in [Zr/Fe] must 
be regarded as well establis hed . The p r oduct ion of Zr in 
AGB stars is discussed by lBusso et al (|2001l ) and more 
recently by iTravaglio et all (|2004fh The yield of Zr via 
the main s-process depends strongly on the size of the 
13 C pocket, as this isotope is the source of the neutrons 
for this channel. The 13 C is itself produced locally within 
the AGB star by proton captures on the abundant iso- 
tope 12 C, but the requirement of having substantial C 
together with protons in a region hot enough for these 
reactions to proceed leads to these reactions occurring 
u nder very constrain ed conditions, as described in detail 
in lBusso et all (|2001l ). The weak s-process, also involved 
in Zr production, ut i lizes 2 2 Ne as the neutron source. 
iKarakas fe Lattariz io (2007) provides an estimate of the 
total yield of elements heavier than the Fe-peak (i.e. be- 
yond the maximum atomic mass included in their reac- 
tion network) in AGB stars. Their results suggest that 
the production of Zr is biased towards lower mass AGB 
stars than those within which the heavy isotopes of Mg 

8 The Ba II lines are too strong to be able to detect the small 
star-to-star variations seen here 
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are made, which could explain why [Zr/Fe] does not ap- 
pear to vary among the M71 giants in our sample. Better 
models to verify the consistency of our results with the 
predicted behavior of Zr and La vs the Mg isotope ratios 
in AGB stars as a function of mass are needed. 

Eu serve the same role for the r-process of neutron 
capture which produces selected isotopes of heavy ele- 
ments beyond the Fe peak. No credible star-to-star vari- 
ations could be detected among the luminous M71 giants 
in [Eu/Fe] based on the strength of the 6645 A line of 
Eu II. 

5. CONCLUSIONS 

With the aid of exquisite spectral resolution and very 
high signal-to-noise spectra combined with a very care- 
ful analysis, we have determined Mg isotopic ratios for 9 
luminous giants in the metal-rich Galactic globular clus- 
ter M71. We have also used these spectra to determine 
precision abundances for several other elements in this 
GC, including the light elements O, Na, Mg, and Al. 

The abundances of Si, Ca, Ti, Ni and Fe do not show 
any star-to-star variations. The total range for the ab- 
solute Fe abundance, log[e(Fe)], among the sample of 9 
giants in M71 is only 0.08 dex. Once the dependence 
on T c ff is removed, all the M71 giants have identical 
[Fe/H], [Si/Fe], [Ca/Fe], [Ti/Fe] and [Ni/Fe] to within 
g = 0.04 dex (10%). This places a strong constraint 
on the uniformity of mixing in the young GC. These ele- 
ments cannot have been produced in anything other than 
the first generation of GC stars. 

We see the expected correlations and anticorrelations 
among the light elements O, Na, Mg, and Al in the M71 
giants. But the amplitude of the star-to-star variations 
among these elements is small, 0.3 dex for [O/Fe], 0.6 dex 
for [Na/Fe J, 0.2 dex for [Al/Fej , and at most 0.1 dex for 
[Mg/Fe]. ICarretta et~al (|2007f) have looked for a link 
between chemical anomalies along the RGB and other 
properties of GCs. In addition to the obvious sugges- 
tion that higher amplitude star-to-star variations should 
be found in higher mass GCs, which with their higher 
binding energies may be better able to retain stellar 
ejecta, they suggest that the high temperature extension 
of the horizontal branch blue tail is longer in GCs with 
higher amplitude Na-0 anticorrelations. The latter is at- 
tributed to a spread in He, and hence may again be tied 
to the ability t o reta in ejected gas, see the discussion in 
ICarretta et all (|2007f ). We suggest that one additional pa- 
rameter is relevant, the Fe-metallicity of the GC. As the 
stellar metallicity becomes higher, the effect of the addi- 
tion of AGB processed material, whose nucleosynthesis 
proceeds among the light elements in a manner which is 
somewhat independent of their initial metallicity, would 
be diluted, and the resulting yield, assumed to be pos- 
itive, reduced. Such an effect is apparent in the ca l- 
culation of AGB yields bv iKarakas fe Lattanz io (2007), 
among others. This would explain why the most promi- 
nent cases of star-to-star variations within GCs are seen 
among the lower metallicity GCs (i.e. M15, M13, etc) 
even though there are a number of quite massive high- 
met al licity GCs such as 47 Tuc for which ICarretta et all 
(|2004h found only modest star-to-star variations among 
the light elements. In this context it is interesting to note 
that the behavior of CH and CN bands in these relatively 
high metallicity GCs tends to be bimodal, while in more 



metal-poor GCs, stars fill the entire range of C and N 
abundances without being so concentrated towards the 
upper and lower extremes of [ C/Fe] (see, e.g. Fig. 11 of 
iCohen. Brilev fe Stetson! 120051) . 

The heavy Mg isotopes among the CN weak giants 
in our M71 sample have very low Mg isotopic ratios 
26 Mg/Mg (~4%) which are consistent with models of 
galactic chemical evolution with no contribution from 
AGB stars. These stars are both CN weak and have 
light element abundances typical of field Galactic halo 
metal-poor stars of similar Fe-meta l licity: we call them 
the "normal" stars. iKroupa fe Boilvl (|2002t ) discuss form- 
ing the Galactic halo field via dissolved GCs; our re- 
sults suggest that this must have involved the "nor- 
mal" GC stars from clusters which dissolved early on 
before their intermediate-mass stars reached the AGB. 
The CN strong M71 giants have higher Mg isotopic ratios 
26 Mg/Mg (~8 %), but these are far below previously pub- 
lished values bv lYong et al.l (|2006f ) and references therein 
for more metal-poor GC giants. This group of stars, in 
addition to a higher fraction of the heavier Mg isotopes, 
shows enhanced Na and Al, accompanied by lower O 
abundances. 

The behavior of the Mg isotopes in the "normal" stars 
is reproduced by m odels of galactic chemical evolution by 
iFenner et all (|2003h without any contribution from AGB 
stars. Their models with the AGB contribution expected 
in the normal course of stellar evolution reproduce the 
behavior of the heavy Mg isotope rich M71 giants. These 
stars must represent a later generation of stars formed 
sufficiently long after the first that the AGB contribu- 
tion to their chemical inventory was included. Alterna- 
tively, the CN-rich and CN-weak stars may have formed 
at the same time, but the CN-weak stars could have been 
formed in an unmixed environment, while the CN-strong 
stars were born from material polluted by AGB stars. 

With the present work, we believe we have demon- 
strated convincingly that the difference between the gen- 
erations of GC stars is consistent in detail with the contri- 
bution of AGB stars. Furthermore, we have ext ended our 
previous results from lMelendez fe Cohenl (|2007l ) based on 
metal-poor halo field dwarfs to still higher metallicity, up 
to the Fe-metallicity of M71. Bas ed on the calculations 
of ste llar yields for AGB stars b y Karak as fe Lattanzio! 
( 2003), updated in IKarakas fe Lattanzio! (|2007l ). stars of 
initial mass 3-6 M Q are required to contribute significant 
amounts of the heavier Mg isotopes. Stellar evolutionary 
tracks establish that stars in this mass range reach the 
upper AGB at an age of ^0.3 Gyr. Assuming uniform 
mixing within the gas in this GC, for a GC as metal- rich 
as is M71, this timescale must reflect the minimum age 
difference between the CN-weak and the CN-rich stars 
in the GC and the maximum age range of the CN-weak 
stars we see today. Furthermore, the enrichment of the 
halo gas from which M71 eventually formed must have 
reached such a high Fe-metallicity without a substantial 
contribution of "polluted" material, whose source we ar- 
gue on the basis of the tight constraints we have placed 
on any variation in abundance among the Fe-peak ele- 
ments in M71 must be AGB stars. 
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TABLE 1 
The List of Atomic Lines Used 



Wavelength 


Exc. Pot. 


log(<?/) 


Species 


(A) 


(eV) 


(dex) 




OoUU.oU 


U.UUU 


— 9. / z 


m Tl 


COCO 77 
DODO. 77 


0.020 


— 10.19 


\r\ Tl 

[O 1J 


6154.22 


2.102 


— 1.547 


Na I 


DlOU. ( 4 


2.104 


— 1.246 


IN a 1 


oolo. ( 1 


5.108 


— 1.945 


TV ,T „ T 

Mg 1 


0319.23 


5.108 


—2.165 


Mg I 


oo9o.ui 


3.143 


— 1.481 


A 1 T 

Al 1 


DD98.DD 


3.143 


— 1.782 


A 1 T 

Al 1 


r /loo no 


5.614 


— 1.69 


Q ■ T 

Ol 1 


CI AO A O 

Dl4Z.4o 


5.619 


1 A 1 

— 1.41 


C3; T 


ooyu.ii 


2.521 


—0.571 


ta 1 


050 i .OD 


2.933 


— l.o / 


T 

ta 1 


oloo.U2 


2.521 


—2.42 


T 

Oa 1 


blbo.4o 


2.521 


— 1.142 


ta 1 


5453.64 


1.443 


— 1.61 


rp- T 
1 1 1 


5648.56 


2.495 


—0.25 


rp- T 
111 


5913.71 


0.021 


— 4.10 


rp- T 
111 


5918.53 


1.066 


— 1.47 


rp- T 
111 


(?nno 7n 

ouy2.7y 


1.887 


— 1.378 


rp- T 
ill 


6273.38 


0.021 


—4.17 


rp- T 
1 1 1 


D7UD.29 


1.502 


—2.78 


rp- T 
1 1 1 


5435.85 


1.986 


—2.60 


INl 1 


5468.10 


3.847 


-1.61 


Nil 


5589.35 


3.898 


-1.14 


Nil 


5846.99 


1.676 


-3.21 


Nil 


6086.28 


4.266 


-0.51 


Nil 


6108.11 


1.676 


-2.44 


Nil 


6130.13 


4.266 


-0.96 


Nil 


6127.45 


0.154 


-1.06 


Zrl 


6134.55 


0.000 


-1.28 


Zrl 


6140.53 


0.519 


-1.41 


Zrl 


6143.20 


0.071 


-1.10 


Zrl 


6390.48 


0.321 


-1.41 


La II 
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TABLE 2 

Stellar parameters, Mg isotopic ratios, C, N, O, C+N, and C+N+O abundances. The scatter (a) of the 25 > 26 Mg/Mg ratios 

OBTAINED BETWEEN THE THREE MgH FEATURES IS GIVEN BETWEEN PARENTHESIS. THE FIRST FOUR STARS ARE CN-STRONG GIANTS AND 

THE OTHER FIVE STARS ARE CN-WEAK GIANTS. 



Star 


T eff ,log g,[Fe/B],v 


m ic ,Vmac , 25 Mg/Mg, 2S Mg/Mg 


[C/Fc] 


[N/Fc] 


[O/Fc] 


A(C) 


A(N) 


A(O) 


A(CN) 


A(CNO) 




(K, dex, dex, km s 


;-\ km s-\ %, %) 


(dex) 


(dex) 


(dex) 


(dex) 


(dex) 


(dex) 


(dex) 


(dex) 


1-45 


3870, 0.60, -0.83, 


1.65, 6.0, 8.4(1.5), 9.2(2.3) 


0.16 


0.75 


0.32 


7.76 


7.75 


8.26 


8.06 


8.47 


I 


4080, 1.20, -0.78, 


1.52, 5.2, 6.8(1.7), 6.9(1.9) 


0.13 


0.60 


0.49 


7.78 


7.65 


8.48 


8.02 


8.61 


1-66 


4150, 1.50, -0.79, 


1.49, 5.6, 6.9(1.9), 7.7(2.5) 


0.12 


0.86 


0.45 


7.76 


7.90 


8.43 


8.14 


8.61 


1-53 


4150, 1.50, -0.78, 


1.52, 5.6, 6.9(0.8), 8.1(1.6) 


0.01 


0.91 


0.40 


7.66 


7.96 


8.39 


8.14 


8.58 


1-46 


3820, 0.45, -0.81, 


1.65, 6.0, 5.9(1.7), 6.0(2.2) 


0.18 


-0.03 


0.39 


7.81 


7.00 


8.36 


7.87 


8.48 


A4 


3870, 0.50, -0.81, 


1.52, 5.5, 3.9(1.7), 5.0(2.1) 


0.17 


0.28 


0.42 


7.79 


7.30 


8.38 


7.91 


8.51 


1-77 


3900, 0.55, -0.80, 


1.57, 5.9, 5.2(1.0), 5.7(1.2) 


0.17 


0.22 


0.40 


7.80 


7.25 


8.37 


7.91 


8.50 


1-64 


4160, 1.40, -0.76, 


1.54, 5.2, 5.9(2.1), 6.0(1.8) 


0.05 


0.38 


0.49 


7.72 


7.45 


8.50 


7.91 


8.60 


1-21 


4270, 1.45, -0.84, 


1.65, 6.0, 5.2(1.9), 3.9(1.6) 


0.19 


0.39 


0.58 


7.78 


7.38 


8.51 


7.93 


8.61 



TABLE 3 

[Na,Mg,Al,Si,Ca,Ti,Ni,Zr,La/Fe] abundance ratios. The FIRST four stars are CN-STRONG giants and the other five STARS ARE 

CN-WEAK GIANTS. 



Star 


[Na/Fc] 


[Mg/Fe] 


[Al/Fc] 


[Si/Fc] 


[Ca/Fe] 


[Ti/Fc] 


[Ni/Fc] 


[Zr/Fe] 


[La/Fc] 




(dex) 


(dex) 


(dex) 


(dex) 


(dex) 


(dex) 


(dex) 


(dex) 


(dex) 


1-45 


0.48 


0.21 


0.36 


0.20 


0.18 


0.29 


-0.02 


0.23 


0.37 


I 


-0.03 


0.15 


0.20 


0.19 


0.16 


0.18 


-0.04 


0.06 


0.29 


1-66 


0.17 


0.14 


0.25 


0.26 


0.16 


0.17 


-0.01 


0.08 


0.36? 


1-53 


0.26 


0.17 


0.28 


0.22 


0.19 


0.22 


0.00 


0.12 


0.43? 


1-46 


0.02 


0.25 


0.26 


0.16 


0.10 


0.27 


-0.06 


0.16 


0.30 


A4 


0.00 


0.24 


0.26 


0.15 


0.18 


0.29 


-0.02 


0.17 


0.29 


1-77 


0.08 


0.19 


0.30 


0.12 


0.22 


0.33 


-0.02 


0.23 


0.28 


1-64 


0.01 


0.18 


0.23 


0.18 


0.20 


0.22 


-0.04 


0.12 


0.29 


1-21 


-0.09 


0.18 


0.18 


0.24 


0.15 


0.12 


-0.04 


-0.05 


0.23 



TABLE 4 

Abundance trends with temperature (A = C + slope * T c ff) for the CN-weak giants 



Abundance 


C 


slope 




(dex) 


(dex K- 1 ) 


[Fe/H] 


-0.7243 


-1.9401E-05 


[C/Fe] 


0.5388 


-9.6626E-05 


[N/Fc] 


-2.4973 


6.8564E-04 


[O/Fe] 


-1.0946 


3.8726E-04 


[Na/Fe] 


0.8767 


-2.1798E-04 


[Mg/Fe] 


0.7654 


-1.3923E-04 


[Al/Fe] 


1.0106 


-1.9097E-04 


[Si/Fe] 


-0.5890 


1.8957E-04 


[Ca/Fe] 


0.0252 


3.6139E-05 


[Ti/Fe] 


1.6833 


-3.5899E-04 


[Ni/Fe] 


-0.0085 


-6.8475E-06 


[Zr/Fe] 


1.9213 


-4.4838E-04 


[La/Fe] 


0.7269 


-1.1193E-04 


A(C+N) 


7.6377 


6.6801E-05 


A(C+N+Q) 


7.3416 


2.9914E-04 




Dh 0.40 



0.00 



• M 71 1-7? (Keck + HIRES, this work) 

■ - M 71 A4 (Subaru + HDS, Yong et al. 2006) 



5134.0 



5136.0 



5138.0 



A (A 



5140.0 



Fig. 1.— Comparison of M71 1-77 (Keck I + HIRES, this work) and M71 A4 (Subaru + HDS, Yong ct al. 2006). Both stars have similar 
stellar parameters (T e ff ~ 3900 K), except that A4 has a lower macroturbulence. Both spectra were taken with a resolving power of f» 10 5 . 




z!.s«\1g/Mg (%) 



28 Mg/Mp = 5.1^0.3 : % 



z"«\lg/'Mg (%) 




z "=M~/\lg (%) 



Fig. 2.— Fits for the 5134.3 A (left), 5138.7 A (center), and 5140.2 A (right) MgH features in the giant M71 1-77. Observed spectra 
are represented as filled circles, and synthetic spectra as solid lines. The calculations were performed for 25 > 26 Mg/Mg ratios of 0-10%. The 
relative variations of the \ 2 fits are shown as a function of the isotopic ratio. 




5131.2 5131.5 5131. B 

MA) 



5 Mg/Mg-\6.9±0.6% 
B Mg/Mg=0.5±O.4% 9 



w.a»Mg/Mg (%) 




5138.6 5138.8 

A(A) 



niH/'Mi4-7;6-0.4% 
a Mg/Mg = 8.^3 = 0.3% 



»Mg/Mg (%) 




5140.0 5M0.2 5140. 

MA) 



Fig. 3.— Fits for the 5134.3 A (left), 5138.7 A (center), and 5140.2 A (right) MgH features in the giant M71 1-53. Observed spectra arc 
represented as filled circles, and synthetic spectra as solid lines. The calculations are shown for 25 > 26 Mg/Mg ratios of 2-12% (left panel) 
and 0-10% (center and right panels). The relative variations of the x 2 fits are shown as a function of the isotopic ratio. 
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Fig. 4. — [Fe/H], [C,N,0/Fe], A(C+N), and A(C+N+0) abundances (dex) as a function of T off . Open and filled circles represent 
CN-weak and CN-strong stars, respectively. The dashed line is a fit to the CN-weak stars, and the dotted lines shows the scatter (<r) 
around the fit. 
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Fig. 5. — 24 > 25 > 26 Mg/Mg ratios and [Na,Mg,Al/Fe] abundance ratios (dex) as a function of T e £f. Open and filled circles represent CN-weak 
and CN-strong stars, respectively. The dashed line is a fit to the CN-weak stars, and the dotted lines shows the scatter (a) around the fit. 
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Fig. 6. — [Si,Ca,Ti,Ni,Zr,La/Fe] abundance ratios (dex) as a function of T e ff . Open and filled circles represent CN-weak and CN-strong 
stars, respectively. The dashed line is a fit to the CN-weak stars, and the dotted lines shows the scatter (a) around the fit. 
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Fig. 7. — 24 Mg/Mg ratios as a function of [Na/Fc], [O/Fc], [Mg/Fe] and [Al/Fc]. Open and filled circles represent CN-weak and CN-strong 
stars, respectively. Trends with T e g have been corrected for the elemental [X/Fe] ratios (O, Na, Mg, Al). 
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Fig. 8. — 25 Mg/Mg ratios as a function of [Na/Fe], [O/Fc], [Mg/Fe] and [Al/Fe]. Open and filled circles represent CN-weak and CN-strong 
stars, respectively. Trends with T e g have been corrected for the elemental [X/Fe] ratios (O, Na, Mg, Al). 
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Fig. 9. — 25 Mg/ 24 Mg ratios as a function of [Na/Fe], [O/Fe], [Mg/Fe] and [Al/Fe]. Open and filled circles represent CN-weak and 
CN-strong stars, respectively. Trends with T c ff have been corrected for the elemental [X/Fe] ratios (O, Na, Mg, Al). 
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Fig. 10. — 26 Mg/Mg ratios as a function of [Na/Fe], [O/Fe], [Mg/Fe] and [Al/Fe]. Open and filled circles represent CN-weak and 
CN-strong stars, respectively. Trends with T e ff have been corrected for the elemental [X/Fe] ratios (O, Na, Mg, Al). 
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Fig. 11. — 26 Mg/ 24 Mg ratios as a function of [Na/Fe], [O/Fe], [Mg/Fe] and [Al/Fe]. Open and filled circles represent CN-weak and 
CN-strong stars, respectively. Trends with T e ff have been corrected for the elemental [X/Fe] ratios (O, Na, Mg, Al). 
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Fig. 12. — Correlations between O, Na, Mg and Al. Open and filled circles represent CN-weak and CN-strong stars, respectively. The 
abundance ratios have been corrected for trends with T e ff . A clear anti-correlation is seen between O and Na, and a well-defined correlation 
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